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ABSTRACT 

We use semi-analytic techniques to study the formation and evolution of brightest 
cluster galaxies (BCGs). We show the extreme hierarchical nature of these objects 
and discuss the limitations of simple ways to capture their evolution. In a model 
where cooling flows are suppressed at late times by AGN activity, the stars of BCGs 
are formed very early (50 per cent at z ~ 5, 80 per cent at z ~ 3) and in many 
small galaxies. The high star formation rates in these high-z progenitors are fuelled by 
rapid cooling, not by merger-triggered starbursts. We find that model BCGs assemble 
surprisingly late: half their final mass is typically locked-up in a single galaxy after 
z ~ 0.5. Because most of the galaxies accreted onto BCGs have little gas content 
and red colours, late mergers do not change the apparent age of BCGs. It is this 
accumulation of a large number of old stellar populations - driven mainly by the 
merging history of the dark matter halo itself - that yields the observed homogeneity 
of BCG properties. In the second part of the paper, we discuss the evolution of BCGs 
to high redshifts, from both observational and theoretical viewpoints. We show that 
our model BCGs are in qualitative agreement with high-z observations. We discuss 
the hierarchical link between high-z BCGs and their local counter-parts. We show 
that high-z BCGs belong to the same population as the massive end of local BCG 
progenitors, although they are not in general the same galaxies. Similarly, high-z BCGs 
end-up as massive galaxies in the local Universe, although only a fraction of them are 
actually BCGs of massive clusters. 

Key words: galaxies: formation - galaxies: evolution - galaxies: elliptical and lentic- 
ular, cD - galaxies: fundamental parameters - galaxies: stellar content 



1 INTRODUCTION 

Brightest cluster galaxies (BCGs) are the very bright galax- 
ies that inhabit the cores of rich galaxy clusters. They are 
among the most luminous and most massive galaxies in 
the Universe at the present epoch. At low redshift, these 
objects exhibit a small dispersion in their aperture lumi- 
nosities (after correction for systematic effects and for de- 
pendences on galaxy structure and environment). As a re- 
sult, many early studies used BCGs as 'standard candles' 
for classical cosmological te sts jSandase|[l973: iGunn fe Okd 
1 19751 : iHoessel fc Schneiderl 1 19851 : iPostman fc Lauerl Il995l ). 
BCGs lie close to the peaks of the X-ray emission 
in concentrated X-ray bri ght clusters (| Jones fc Formanl 
1 1984 iRhee fc Latouil Il99lh and their rest-frame veloc- 
ities have small offsets from those of their host clus- 
ters (iQuintana fc Lawrielll982l ; IZabludoff. Huchra fc Gellerl 
1990). They do not appear to be dra wn from the same lumi- 
nosity function as cluster ellipticals (|Tremaine fc Richstond 
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1 19771: iDressien 1978) or a s bright galaxies in general 
( Bernstein fc Bhavsarl 1200 ll ). They also exhibit different 
luminosity profiles than typic al cluster elliptical galaxies 
(|Oemlerll 19761 : ISchombertll 19861 ) . All these observations sug- 
gest that the evolution of the BCGs might be appreciably 
distinct from galaxy evolution in general. 

Early theoretical studies dis cussed the role of X- 

ray d riven c o oling flows (Silk] Il976l : ICowie fc Binnevl 
1 19771 : iFabianl 1 1994ft and cannibal i sm d u e to dynam- 
ical friction ifostriker fc Tremainel Il975l: IWhitel Il976l : 
iMalumuth fc Richstond Il984l : lMerrittJll985h in the forma- 
tion of these special objects. This early work was not how- 
ever successful due to the use of a simplified cluster model. 
In the now standard cold dark matter model of structure for- 
mation, we understand that local massive haloes (clusters) 
assembled rather late, through the merging of smaller sys- 
tems. In this perspective, cooling flows are indeed the main 
fuel for galaxy mass-growth at high redshifts, in dense and 
lower mass haloes. This source is removed only at low red- 
shifts and in group or cluster environments, possibly due to 
AGN feedback. In addition, galaxy-galaxy mergers are most 
efficient within small haloes with low velocity dispersions. 
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They are indeed driven by dynamical friction, but it is the 
accretion rate of galaxies into the proto-cluster, along with 
cluster growth itself, that regulates and sets the conditions 
for galaxy merging. 

Since these early studies, little theoretical work has fo- 
cused on the formation and evolution of BCGs, and it was 
only in 1998 that t he subject was revisited in a CDM hierar- 
chica l framework (Dubinski 1998; Aragon-Salamanca et all 
1 19981 ). iDubinskil l|l998h used a N-body simulation of a 
cluster of galaxies and showed that merging naturally 
produces a massive central galaxy with surface bright- 
ness and velocity dispersion that resemb le those of BCGs. 
lAragon-Salamanca. Baugh fe Kauffmannl |l998) analysed 
the K-band Hubble diagram of BCGs up to redshift ~ 1 
and interpreted the apparent lack of passive evolution as ev- 
idence for recent mass accretion by the BCGs. Using a sim- 
ple parametrisation, the authors estimated a mass growth 
by a factor 4 — 5, a result that they found to be consis- 
tent with predictions fr om semi-analytic models of galaxy 
formation. Later work (Col lins fc Mannlll998l ; iBurke et al.l 
2000; iBrough et al.l 120021 ) questioned th is agreement show- 
ing that it might reflect clu ster selection. Burke et al.l (2000) 
showed that the results of I Aragon-Salamanca et al.l ( 1998 ) 
were based on a sample of low X-ray luminosity (low-Lx) 
clusters. For high-Lx clusters, which better match the clus- 
ter mass selection used in the semi-analytic models, they 
found substantially lower rates of mass accretio n, at most 
a fact or of ~ 2 since z ~ 1. In recent work, iGao et all 
(|2004l ) used numerical simulations to study the formation 
of the inner cores of massive dark matter haloes. Their re- 
sults suggest that typical central galaxies have undergone a 
significant number of mergers since z ~ 1. Observations do 
prov ide some examples of central galaxies in the act of merg- 
ing l|Nipoti et al.l [2003h . although not many. It is therefore 
still a matter of debate whether the large number of merg- 
ers predicted in a hierarchical scenario is in agreement with 
observational data. 

In this paper we study the formation and evolution of 
BCGs using a combination of N-body simulation and semi- 
analytic techniques. In these models, central galaxies are 
bound to be 'special': gas that is shock heated to the virial 
temperature of dark matter haloes cools radiatively and con- 
denses only onto the central galaxy. An additional physical 
process, such as feedback from a central AGN, must be in- 
voked to prevent excessive growth from such cooling flows. 
The haloes which BCGs inhabit originate from the gravita- 
tional collapse of the highest (and rarest) peaks of primor- 
dial density fluctuations. In these overdense regions, haloes 
collapse earlier and merge more rapidly in comparison to re- 
gions of the Universe with average density. Taking advantage 
of the largest simulation of the cosmi c structure carried ou t 
so far - the Millennium Simulation l|Springel et alj|2005h - 
we select a large number of massive haloes at z = and 
study in detail how and when the stars that compose their 
central galaxies formed, and how the assembly of these ob- 
jects relate to the assembly of the halo itself. In the second 
part of the paper, we investigate if the low rate of mass ac- 
cretion inferred from observations is in agreement with our 
model predictions. 

We describe the simulation and the semi-analytic model 
used for this study in Sees. [2] and [3] In Sec. [4] we present a 
case-study in order to define the terminology that we will use 



later on. In Sec.[5]we present statistical results for 125 BCGs 
selected from the simulation at z — 0. In Sec. [6] we extend 
that selection to higher z and we compare the luminosity 
evolution measured using high-Lx clusters to that predicted 
by our model. Finally, we discuss our results and give our 
conclusions in Sec. [7] 



2 SIMULATION AND HALO MERGER TREES 

In this study, we make use of the Millennium SimulatioiQ 
recently carr ied out by the Virgo C onsortiunjfl and described 
in detail in ISpringel et all (|2005D . The simulation follows 
N = 2160 3 particles of mass 8.6 x 10 s h~ 1 M Q within a co- 
moving box of size 500 /i _1 Mpc on a side. The cosmological 
model is a ACDM model with fi m = 0.25, fi b = 0.045, 
h — 0.73, = 0.75, n = 1, and ag = 0.9, where the Hub- 
ble constant is parameterised as Ho = 100 h kms _1 Mpc _1 . 
These cosmological parameters are consistent with recent 
determinations from the co mbined analysis of th e 2dFGRS 
and first year WMAP data (|Sanchez et all 120061 ). Given its 
high resolution and large volume, the Millennium Simulation 
allows us to follow in enough details the formation history 
of a representative sample of rare massive clusters. 

During the simulation, 64 snapshots were saved, 
together with group catalogues and their embedded 
substructures, ident ified using the algorithm SUBFIND 
j|Springel et alj|200ll ). As explained in lSpringel et all j|2005l ). 
dark matter haloes are identified using a standard friends- 
of-friends (FOF) algorithm with a linking length of 0.2 in 
units of the mean particle separation. The FOF group is 
then decomposed into a set of disjoint substructures, each 
of which is identified by SUBFIND as a locally overdense re- 
gion in the density field of the background halo. Only the 
identified substructures which retain after a gravitational 
unbinding procedure at least 20 bound particles are consid- 
ered to be genuine substructures. We note that SUBFIND 
classifies all the particles inside a FOF group either as be- 
longing to a bound substructure or as being unbound. The 
self-bound part of the FOF group itself will then also appear 
in the substructure list and represents what we will refer to 
below as the main halo. This particular subhalo typically 
contains 90 per cen t of the mass of the FOF group (e.g. 
ISpringel et aLll200ll ). The group catalogues were then used 
to construct detailed merging history trees of all gravitation- 
ally self-bound dark matter structures. These merger trees 
form the basic input needed by the semi-analytic model used 
here. 

In the present paper, we use an improved scheme for 
tracking halo c entral galaxies, which we briefly explain in 
the following. In lSpringel et al.1 l|2005l) . the "first progenitor" 
of a given halo was simply defined as the most massive of 
its progenitors. This pointer was meant to efficiently track 
the main branch of a merger tree, i.e. the branch containing 
the central galaxy of a FOF group. However, this selection 
is sometimes ambiguous, for example when there are two 
subhalos of similar mass such that their rank order in size is 
subject to noise. In order to avoid such occasional failures, 
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we have redefined the most-massive progenitor branch by 
using a different criterion to determine the "first progenitor" 
designation. Briefly, let Ni denote the number of self-bound 
particles of halo i. Then we define a quantity Mi as 



Mi = Ni+ max(M il , Mi, 



(1) 



where the maximum is taken over all the progenitors ij of 
halo i. This definition is applied recursively to each halo, 
over all timesteps, and sets Mi to the mass of the main 
branch rooted in halo i. The "first progenitor" is then se- 
lected as the progenitor with the largest value of Mi. This 
definition therefore selects, for any halo, the branch that ac- 
counts for most of the mass of the final system for the longest 
period. We argue that this provides a reasonable and very 
robust definition of the "most ma ssive progenito r histo ry" . 
Note that the scheme adopted in ISpringel et all (|2005l ) for 
selecting the main progenitor branch corresponds to replac- 
ing equation (fTJ with Mi — Ni. 



3 SEMI-ANALYTICS 

The semi-analytic model we use in the pr esent paper is 
a slig ht ly modified v e rsion of that used in Springel et al.l 
(2005), ICroton et all |2006t ). and |Pe Lucia et al.1 (|2006). 



This model builds on the m ethodology int r oduce d by 
Kauffmann fc Haehneltl d200Ch . ISpringel etail <|200ll ) and 
De Lucia. Kauffmann fc White! l|2004l ), and includes a 



model for the suppre ssion of cooling flow s by "radio mode" 
feedback from AGN (|Croton et alj|2006h . The reader is re- 
ferred to these papers for a detailed description of the phys- 
ical processes included in the model. In this section, we 
simply summarize our additions and changes relative to 
ICroton et all (|200tf) (Sees. E2] and G3J, and remind the 
reader with our treatment of some of the processes that are 
relevant for the present work (Sec. I3.1|l . The prescriptions 
and parameter values described below are the ones that have 
been used to generate our standard mo del accessible online 
|Lemson fc the Virgo Consortium! 120061 ). 



3.1 Galaxy Mergers 

We treat galaxy m ergers as in ISpringel et al.l (|200lh and 
ICroton et al.l l|2006l ). Substructures allow us to follow prop- 
erly the motion of the galaxies sitting at their centres 
until tidal truncation and stripping disrupt the subha- 
los at the reso lution limit of the s imulation (here 1.7 x 



10 w M^h- 1 ) jGhigna et~ai1 |200C|; jKravtsov et ail |2004 



|Pe Lucia et al.ll2004l ; iGao et al.|[2004) . When this happens, 
we estimate a survival time for the galaxies using their cur- 
rent orbit and the c lassic al dynamical friction formula of 
iBinnev fc Tremaind l| 19871 ). After this time, the galaxy is 
assumed to merge onto the central galaxy of its own halo. 
This is usually the main halo of the FOF group but can 
be a proper substructure within it. In this paper, we have 
opted to increase the merging time estimates by a factor 
of 2 compared to t he dynamical friction formula used by 
ICroton et al.l ([2006) . This change is not essential for the 
success of the model but it slightly improves the fit of the 
bright end of the luminosity function. The adjustment is 
somewhat ad-hoc but seems justified given that the terms 
that enter the estimation of the dynamical friction times 



carry large uncertainties (most notably the Coulomb loga- 
rithm and the orbital distribution) , and that the correct pre- 
factor has not been calibrated yet with detailed numerical 
studies at the resolution we have here. In addition, there are 
some dire ct indications that th ere are problems with these 
formulas: ISpringel et al.1 (|200ll ) noted that when haloes of 
comparable mass merge, the inferred merging times are typ- 
ically shorter than the time measured using high resolution 
numerical simulations. 

A galaxy merger is accompained by a starburst mod- 
elled usi ng the "collisional star burst" prescription intro- 
duced by ISomerville et al.l j|200ll ). Say two galaxies Qi and 
Qi of masses mi > m.2 merge together. We assume that all 
the gas from Q\ and Q2 is gathered in the disk component of 
the remnant galaxy Q, while all the stars of Q2 are added to 
the bulge stars of Q\ into the buge component of Q. A frac- 
tion of the gas in Q is then converted instantaneously into 
stars, the fraction depending on the baryonic mass ratio of 
the two merging galaxies: 



0.56 x 



Vmi/ 



x m ga 



The numerical parameters in the eq uation above prov ide a 
good fit to the numerical results of ICox et al.l (|2004l ). We 
note that the above formulation implies that in the case of 
an equal mass merger (mi ~ m.2), about 40% of the gas is 
turned into star^f], while in the case where 7712 << mi, no 
star formation is triggered by the interaction. In case of a 
major merger, which occurs when 7x12/ m\ > 0.3, we further 
assume that the final disc is completely disrupted and all 
the stars of galaxy Q are put in a bulge component. 

The formation of new stars is accompanied by the ejec- 
tion of a fraction of the remaining cold gas, which is mod- 
elled u sing the same prescription as adopted in lCroton et al.l 
(2006, see their Sec. 3.6). Note th at this feedback mode l 
is different than that employed in |Pe Lucia et al.l (|2006l ) 
and results in more efficient ejections especially for massive 
galaxies. For the purposes of this work, these two schemes do 
not produce significantly different results, except for Figs. 1131 
and 1141 which we discuss in Sec. 16.11 

3.2 Stellar Populations 

The photometric properties of our model galaxies are com- 
puted employing th e stella r population synthesis model from 
IBruzual fc Chariot! (12003!) and using t he method describe d 
De Lucia et al.l (|2004T l. Contrary to lCroton et al.l (l2006h, 



we ad opt here the initial mass function (IMF) from I Chabrier! 
(2003), together with the Padova 1994 evolutionary tracks. 
The Chabrier IMF can be represented by a power-law with 
index 1.3 from 1 to 100 Mq, and a lognormal distribution 
from 0.1 to 1 Mo, centered on 0.08 Mq and with dispersion 
0.69 ("see IBruzual fc Charlo"tll2003l ). We note that the spec- 
tral properties obtained using this IMF are very similar to 
those obtained using the Kroupa (2001) IMF. In particular, 
a Chabrier IMF yields a larger fraction of massive stars per 
unit mass than Salpeter IMF, and hence a mass-to-light ra- 
tio typically about 1.8 times lower in the V band (see Fig. 

3 The fraction of gas that is turned into stars during a major 
merger saturat es at 40% because of the competition of feedback, 
as explained in lCroton et al. I ||2006|) . 
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4 of iBruzual fe Charlotl lioblh . The Chabrier IMF is physi- 
cally motivated and provides a better fit to coun ts of low- 
mass stars and brown dwarfs in the Galactic disc l)Chabrierl 
I2003T I. 

Because of the substantial differences between the 
Salpeter and Chabrier IMFs, we had to re-adjust some of the 
model parameters in order to obtain a good agreements with 
observati onal results in the lo cal Universe. With respect to 
Table 1 in lCroton et alJ (|2006h . we had to increase kagn (the 
quiescent hot gas BH accretion rate) from 6 x 10 -6 Moyr -1 
to 7.5 x 10 -6 MQyr -1 and decrease osf (the star formation 
efficiency) from 0.07 to 0.03. Finally, for consistency with the 
use of a Chabrier IMF, we adopt an instantaneous recycled 
fraction R = 0.43. 

3.3 Dust attenuation 

In this work we use a new parametrisation for 
dust attenuation, which combines those developed in 
iDevriendt. Guiderdoni fc Sad at (1999) - for a homogeneous 
inter-s tellar medium (ISM) component - and lCharlot fc FaDI 
(2000) - for molecular clouds around newly formed stars. 

We assume that the mean perpendicular optical depth 
of a galactic disk at wavelength A is: 



T\ = 



2.1 10 21 atoms cm" 



where the mean H column density is given by: 

Af co ld 



1.4m p 7r(a 7-1/2) 



-atoms cm 



In the previous equation, ri/2 is the half-mass radius of the 
disk and a = 1.68 is such that the column density repre- 
sents the mass-weighted average column densi ty of the disk, 
which is as sumed to be exponential. As in iHatton et al.l 
|2003h - see lGuiderdoni fc Rocca-Volmerangel (|l987T ) for de- 
tails - the extinction curve depends on the gas metallic- 
ity Z g and is based on an interpolation between the Solar 
neighbourhood and the Large and Small Magellanic Clouds 
(s = 1.35 for A < 2000A and s = 1.6 for A > 20001). 
T he adopted extinction curve for so lar metallicity is that 
of iMathis. Mezger fc Panagial (1983). Finaly, we assign a 
random inclination to each galaxy and apply the dust cor- 
rection to its disc co mponent using a 'slab' geometry (see 
IDevriendt et al.lll999l ). 

In addition to this extinction from a diffuse ISM compo- 
nent, we have also implemented a simple model to take into 
account the att enuation of young sta rs within their birth 
clouds, based on ICharlot fc FaU (|2000h . Stars younger than 
the finite lifetime of stellar birth clouds (which we assume to 
be equal to 10 7 yr) are subject to a differential attenuation 
with mean perpendicular optical depth: 



and 



V 5500^1/ 



BC z „ I 1 , 
TV = T\ X [ 1 



where n is drawn randomly from a Gaussian distribution 
with centre 0.3 a nd width 0.2, truncated at 0.1 and 1 (see 
iKong et alj|2004| ). 



4 A CASE STUDY 

In the framework of a hierarchical scenario, the history of 
a galaxy is fully described by its complete merger tree. 
Whereas in the monolithic approximation the history of a 
galaxy can be described by a set a functions of time, hierar- 
chical histories are difficult to summarise in a simple form, 
because even the identity of a galaxy is sometimes ill-defined. 
In hierarchical models, "a galaxy" should not be viewed as 
a single object evolving in time, but rather as the ensem- 
ble of its progenitors at each given time. In this section, wc 
focus on a single model BCG and describe the formation 
and assembly history of its stars in detail. This allows us to 
define and illustrate the behaviour of quantities that cap- 
ture different aspects of BCG evolution. We will use these 
quantities later (in Sec. [SJ to describe BCG formation in a 
statistical fashion. 



4.1 Trees and branches — definitions 

Fig. [T] shows the full merger tree of the central galaxy of a 
dark matter halo of mass M200 = 8.9 x 10 14 Mq at Z = 0. 
The BCG itself lies at the top of the plot (at z = 0), and all 
its progenitors (and their histories) are plotted downwards 
going back in time recursively. Galaxies with stellar mass 
larger (resp. smaller) than 10 10 /i _1 Mg are shown as sym- 
bols (resp. lines), and are colour-coded as a function of their 
rest-frame B - V colour. The left-most branch in Fig. [T]rep- 
resents what we will hereafter refer to as the main branch. 
This particular branch is obtained by connecting the galaxy 
at each time-step to the progenitor with the largest stellar 
mass (the main progenitor) at the immediately preceding 
time-step. It is tempting to consider this branch as the BCG 
itself, and the objects that merge onto it as other progenitors 
of the BCG. A clear distinction between "main" and "other" 
progenitors is indeed often quite appropriate - typically for 
late type galaxies - and it is satisfactory as long as galaxies 
merging onto the main branch have stellar masses that are 
much smaller than the current mass of the main progenitor. 
In this case, the evolution of the BCG can be characterised 
by a series of accretion events (of much smaller objects) that 
do not introduce major changes in the stellar population or 
identity of the BCG itself. 

This simple summary of a merger tree is however not 
sufficien t to describe historie s as complex as those of BCGs 
(see also iNeistein et al"1l2006l ). Fig. [1] shows that in our case- 
study, the main branch captures the evolution of the BCG 
itself for the last ~ 8 Gyr. Before that time, the main pro- 
genitor of the BCG is only marginally more massive than 
progenitors in other branches. At this point, choosing one 
branch or another becomes arbitrary and a single branch is 
certainly a poor proxy for describing the evolution of the 
BCG and of its stellar population. In the following, we will 
simply refer to the main progenitor of the BCG at any given 
time as the main branch without implying that it necessarily 
contains most of the stars. 

Galaxies that merge onto the main branch must first be 
accreted onto the same halo, and it is therefore interesting 
to establish the connection between the galaxy and the halo 
merger trees. Fig. shows the full tree of the FOF-group 
containing our case-study BCG. The branch highlighted in 
green is the branch containing the main branch of the BCG. 
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Figure 1. BCG merger tree. Symbols are colour-coded as a function of B - V colour and their area scales with the stellar mass. Only 
progenitors more massive than 10 10 Mq ft -1 are shown with symbols. Circles are used for galaxies that reside in the FOF group inhabited 
by the main branch. Triangles show galaxies that have not yet joined this FOF group. 




Figure 2. Merger tree of the FOF group in which the BCG sits at redshift zero. Only the trees of subhalos with more than 500 particles 
at z = are shown. Their progenitors are shown down to a 100 particle limit. Symbol coding is the same as in Fig. [T] The left-most tree 
is that of the main subhalo of the FOF, while the trees on the right correspond to other substructures identified in the FOF group at 
z = 0. In green, we mark the subhalo that contains the main branch of the BCG. 
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The right-most branches are merger trees of secondary sub- 
structures (only those with more than 500 particles) present 
in the FOF-group. These substructures have not yet dis- 
solved into the main halo, and their galaxies can thus not 
contribute to the merger tree of the BCG. In Figs. [T] and 
[5] circles mark objects (galaxies or haloes) that belong to 
the same FOF group as the main branch of the BCG, while 
triangles mark objects that have not yet joined the FOF 
group. Typically, when a halo is accreted onto a bigger sys- 
tem (i.e. joins the same FOF group), it loses mass efficiently 



0.2 



redshift 
0.5 1 



3 5 11 



due to tidal stripping (|Ghigna et~af 2000l ; iKravtsov et al.l 
l2004lDe Lucia et alJl2004lGao et al.l2004h . A nice example 
of this process is shown by the halo branch located roughly 
at the centre of Fig. [2] It is only when the subhalo dissolves 
that its galaxies become part of the main halo of the FOF 
group and are then allowed to merge with the central galaxy 
on a dynamical friction timescale. 

Given the complexity of the merger history shown in 
Fig. [1] it is helpful to define several times that mark impor- 
tant phases in the evolution of a BCG. We call identity time 
(tid) the cosmic time when the BCG acquires its final iden- 
tity. We define tid as the time when the last major merger 
on the main branch occured, i.e. when the most massive 
galaxy which merges on the main branch is more massive 
than a third of the mass of the main progenitor. Before tid, 
the BCG does not exist as one single object, but as sev- 
eral progenitors of comparable masses. Our definition of the 
identity time can be extended to account for multiple simul- 
taneous mergers. We thus define the extended identity time 
(tid) as the latest (cosmic) time when the sum of the masses 
of the progenitors merging on the main branch was larger 
than a third of the mass of the main progenitor. By defini- 
tion, tid < tid, although they are equal in most cases (see 
Fig. 1 in Sec. . 

Early theoretical work discussed the difference between 
formation and asse mbly times for elliptical galaxies in a hi- 
erarc hical context (|Baugh. Cole fc Frenkl [19961 : iKauffmannl 
1996), although this differen ce has been quantified only re- 
cently (|De Lucia et al.l 1200 6). Following this work, we call 
assembly time (t a ) the time when the main progenitor con- 
tains half the final stellar mass of the BCG. As discussed 
above, although the main branch may capture the identity 
of a galaxy for some time, it is not suited to describe the 
complete evolution of the BCG. The stellar population of 
the BCG, in particular, can be fully described only by taking 
into consideration the whole tree, because a large fraction 
of stars actually form in secondary branches. It is therefore 
useful to define a more classical formation time (tf) as the 
time when the total mass of stars formed reaches half the 
final mass of the BCG. By "total mass" we mean, at each 
cosmic time, the sum of the stellar masses of all the progen- 
itors present at that given time, i.e. the projection of Fig. [1] 
onto the vertical time axis. 



4.2 Mass build-up of the BCG 

Fig. Q] shows that the stars that end up in the BCG today, 
start forming at very high redshifts. Rapid cooling in the 
early phases of the cluster collapse lead to the formation 
of a massive central galaxy of stellar mass ~ 
about 11 Gyr ago (z ~ 2.5). A number of accretions of mas- 
sive satellites increases its mass to slightly less than half its 
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lookback time (Gyr) 



14 



Figure 3. Formation and assembly histories of our case-study 
BCG. The black (resp. green) line shows the stellar mass in the 
main branch (resp. the mass summed over all progenitors) of our 
case-study BCG versus time. The horizontal dashed line corre- 
sponds to half the present mass of the BCG. The vertical lines 
mark the identity and extended identity times (see text for de- 
tails). 



present value at redshift about 1 and then the BCG con- 
tinues growing by accretion of satellites. Fig. Q] also clearly 
shows that the satellites accreted below redshift ~ 1 are red 
and much less massive than the main branch. 

In Fig. [3]we show the 'formation' and 'assembly' histo- 
ries of the stars that end up in the BCG of our case-study 
cluster. The black line in Fig. [3] shows the stellar mass of 
the main branch. The green line shows the sum of the stellar 
masses in all progenitors at each time. The horizontal dashed 
line corresponds to half the present mass of the BCG. This 
figure clearly shows that while the main progenitor reaches 
half its final mass only at redshift z a ~ 0.5, half of the stars 
are already formed at redshift zj ~ 5. The identity time and 
the extended identity time are equal in this particular case 
and are indicated by the vertical line at z ~ 1.2. 

It is now interesting to investigate in more detail how 
the BCG gained its mass and how the stars that end up in 
it were formed in the first place. Was star formation mainly 
triggered by starbursts in an early and dense environment? 
Or was star formation mostly quiescent, fuelled by rapid 
cooling? 

Fig. [4] shows the mass in the main branch as a func- 
tion of time (solid line), and the fraction of it that was 
gained through accretion of galaxies (dash-3-dot curve). 
The dashed (resp. dot-dashed) curve indicates the contri- 
bution to accretion by galaxies with stellar masses lower 
(resp. larger) than 10 10 /i _1 Mq (the same cut as in Fig. [1]). 
Fig. [4] shows that most of the stars of this BCG were ac- 
tually not formed in the main branch, but were instead ac- 
creted steadily over time. It is also interesting to note that, 
although the BCG tree has many branches, only a small 
fraction of these contribute significantly to the build-up of 
its mass: ~ 70 per cent of the mass comes from accretion 
of 12 galaxies more massive than 10 10 h~ 1 M.Q. Note also 
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Figure 4. Accretion history of our case-study BCG : fraction 
of the stellar mass that was brought by mergers onto the main 
branch. The solid line is the mass of the main branch. The 3- 
dot-dashed (resp. dot-dashed, dashed) line is the mass accreted 
by all (resp. more massive than 10 10 /i _1 Mq, less massive than 
10 10 ft- 1 M ) galaxies merging on the main branch. 



that the main progenitor becomes rapidly massive enough 
that most galaxies falling onto the main branch give rise 
only to minor merger events. Fig. [4] shows that since z ~ 4 
the mass growth of the main branch is only due to accre- 
tion: AGN activity has switched off cooling and hence star 
formation in the main branch. The stars that compose the 
BCG today formed in separate entities, mostly as a result 
of quiescent star formation: less than 10 per cent of the 
stars were forme d during starbursts (in agreeme nt with pre- 
vious results by lAragon-Salamanca et al.lll998h . Note that 
this quiescent-to-bursty star formation ratio depends on the 
st ar formation prescription that we use for di s cs. Th e results 
of ISomerville et al.1 i|200lj ) and iBaugh et alj (|2005l ) suggest 
that our model predicts an upper boundary of this ratio. 

In a hierarchical framework, the formation of a BCG 
naturally relates to the formation of the cluster itself. Fig. [S] 
shows how the stellar mass in the main branch (in orange) 
compares to the evolution of the dark matter (solid black 
line) and stellar content (dashed black line) of the halo in 
which the main progenitor of the BCG resides at each time. 
The red lines show instead the evolution in dark matter 
(solid) and stellar content (dashed) of all progenitors of the 
halo in which the galaxy sits today. Dark matter masses are 
normalised to the dark matter mass of the FOF group at 
redshift zero, while stellar masses are normalised to the to- 
tal stellar content of the FOF group at redshift zero. The 
figure shows that the assembly history of the BCG follows 
the evolution of the dark matter and stellar content of its 
own halo with a delay due to dynamical friction. The par- 
ent halo of the BCG represents only 90 per cent of the total 
mass of its parent FOF group, and its total stellar content 
accounts for about 35 per cent of the stellar content of the 
FOF group. Note that although dark matter substructures 
only represent ~ 10 per cent of the mass of the FOF group, 
they contain about 65 per cent of its stellar content. Inter- 
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Figure 5. Evolution of the stellar content (dashed lines) and dark 
matter content (solid lines) of the subhalo in which our case- 
study BCG sits at z = 0. Black and red lines are obtained by 
following the main progenitor only and all progenitors at any 
given redshift respectively. Dark matter masses are normalised to 
the dark matter mass of the FOF group at redshift zero, while 
stellar masses are normalised to the total stellar content of the 
FOF group at redshift zero. The orange line shows the stellar 
mass in the main branch. 



estingly, the total stellar content of all the progenitors of the 
halo increases very rapidly and stays constant after redshift 
~ 1. This indicates that many substructures accreted below 
redshift ~ 1 survive as independent entities to the present 
day so that their galaxies do not contribute to the stellar 
content of the halo. 



5 STATISTICAL RESULTS 

In the previous section we have analysed in detail the for- 
mation and assembly history of one BCG chosen as a case 
study, and we have related these histories to the history of 
the parent halo itself. In this section, we extend this anal- 
ysis to a statistical sample of BCGs selected from our sim- 
ulation and study the dispersion in the quantities defined 
above. This is of particular interest because observed BCGs 
exhibit a remarkably small dispersion in their luminosities 
(see Sec. D and stellar popu l ation at low redshift (e.g . 
iMcNamara fc O'Connell Il989l : |j ames fc Mobasher] [2000). 
This is commonly interpreted as a clear indication that these 
objects had similar formation histories. In order to address 
this point, we have selected the 125 haloes more massive 
than 7 x 10 14 M Q at z = from the Millennium Sim- 
ulation and analysed their formation and assembly histo- 
ries, as in the previous section. The BCGs of these objects 
have mean absolute magnitude Mk = —26.58 and a disper- 
sion of 0.20 mag. These values appear to be in nice agree- 
ment with observation al results (see for example Fig. 10 in 
ICollins fc Manrilll998h . This agreement represents a major 
success of the underlying galaxy formation model. 

In Fig. [6] we show the fraction of clusters with identity 
(green), extended identity (red), assembly (blue), and for- 
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Figure 6. Fraction of clusters with identity (green), extended 
identity (red), assembly (blue) and formation (orange) lookback 
times larger than that shown on the x-axis. 



mation (orange) lookback times larger than that shown on 
the x-axis. The green line in Fig.[6]shows that for only about 
40 per cent of our BCGs was the last major merger before 
redshift 1, and for ~ 65 per cent of them it occurred before 
redshift ~ 0.5. The numbers are slightly lower if one con- 
siders, instead of the time corresponding to the last major 
merger, the time when the sum of the masses of the accreted 
objects is equal to one third the mass in the main branch at 
the time of accretion (the extended identity time). The iden- 
tity of the main progenitor of a redshift zero BCG is thus 
typically ill-defined before z ~ 0.7. The assembly times mea- 
sured for our BCGs are low : half of them assembled after 
z ~ 0.5, and only 10 per cent of our BCGs assembled before 
z ~ 1. In contrast, for almost 100 per cent of our BCGs, half 
of the stars were already formed at redshift larger than 4, 
which indicates that all our BCGs have uniformly old stellar 
populations. 

Thick lines in Fig. [7] show the median assembly (blue) 
and formation (green) histories of our sample of redshift 
zero BCGs (as in Fig. [3)l . The dashed regions in this figure 
indicate the 15th and 85th percentiles of the distributions. 
Interestingly, there is a very small scatter in the formation 
histories of the BCGs: for essentially all the objects in our 
sample, 50 per cent of the stars are already formed at red- 
shift 5 (as noticed from Fig. [6} ■ The assembly histories ex- 
hibit a much larger scatter with the fraction of mass in the 
main progenitor varying between 15 and 40 per cent at red- 
shift 1, and between 40 and 70 per cent at redshift 0.5. In 
Sec. [4] we have seen that the mass of the BCG mainly in- 
creases through accretion of satellites, Fig.[7]shows that this 
is the case for the whole sample. The difference between for- 
mation time and as sembly time that we fi nd here is higher 
than that found by |Pe Lucia et all l|2006h for ellipticals, as 
expected from the "extreme" nature of BCGs. We note also 
that the we find a median mass growth of a factor ~ 3 from 
z ~ 1 to 2 = 0. This is in apparent contradiction with ob- 
servational data mentioned in Sec.[T] We will show that this 
is probably only an 'apparent' contradiction in Sec. 16.11 
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Figure 7. Assembly (blue) and formation (green) histories of our 
sample of BCGs selected at redshift (as in Fig. [3). Thick lines 
show the median of the distributions, while the dashed regions 
show the 15th to 85th percentile range. 
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Figure 8. Evolution of the dark matter mass of the FOF group 
containing the main progenitor of the BCG at each time (blue), 
of the stellar content of this FOF group (green), and of the mass 
of the main progenitor of the BCG (orange). Thick lines show 
the median and shaded regions show the 15th to 85th percentile 
range. Note that we have multiplied the stellar mass of the FOF 
group by a factor of 50 and the mass of the main progenitor of 
the BCG by a factor of 100. 



Fig. [5] shows how the evolution of the mass in the main 
progenitor (orange) compares to the evolution in the stellar 
(green) and dark matter (blue) content of the FOF group 
in which the main progenitor resides at each lookback time. 
Thick lines show the median and shaded regions show the 
15th and the 85th percentiles of the distributions. Note that 
we have multiplied the stellar mass of the FOF group by a 
factor of 50 and the mass of the main progenitor of the BCG 
by a factor 100 for display purposes. The figure shows that 
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Figure 9. Distribution of objects accreted onto the main branch 
as a function of their gas fraction (panel a), star formation rate 
(panel b), stellar metallicity (panel c), and B-V colour (panel d). 
Green histograms refer to all progenitors, independently of their 
accretion redshifts. Red histograms refer to objects accreted at 
redshift larger than 3. 
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Figure 10. Morphological mix of objects accreted onto the main 
branch in different redshift bins. The red indicates early-type mor- 
phologies, the green is for SO and the blue is for late types (see 
text for details). For each redshift bin, the left-most histogram 
refers to a normalisation in mass, while the right-most histogram 
refers to a normalisation in number. 



the stellar content of the FOF group in which the BCG sits, 
represents about 1.7 per cent of its dark matter mass, at 
all redshifts. In turn, the mass in the main branch is ~ 100 
per cent of the stellar content of the FOF group at high z, 
and then decreases to reach a plateau at z ~ 1, with a mass 
fraction of ~ 6 per cent. 

The example in Fig. [1] shows that most of the progeni- 
tors accreted onto the main branch, are already red at the 
time of accretion. We now want to quantify the distribution 
of the objects accreted onto the main branch in terms of dif- 
ferent physical properties. In Fig. |5]we show the distribution 
of these progenitors as a function of their gas fraction (panel 
a), star formation rate (panel b), stellar metallicity (panel 
c), and B-V colour (panel d). Note that we are using here 
only progenitors more massive than 10 10 /i _1 M0, as these 
are the ones that contribute most significantly to the mass 
build-up of the BCG (see Fig. [3} . The green histograms in 
all panels refer to all the objects accreted, independently of 
their accretion redshifts. The top left panel indicates that 
the great majority of these progenitors have very low gas 
fractions. Only objects accreted at very high redshift have 
some residual gas, as indicated by the red histograms; these 
refer to progenitors accreted at redshift larger than 3. As a 
consequence, only objects accreted at very high redshift have 
some residual star formation (panel b) and relatively blue 
colours (panel d). Interestingly, the metallicity distribution 
of progenitors does not show any significant evolution with 
redshift. It peaks at Z ~ 0.5Zq and is quite broad. 

Fig. [10] shows the morphological mix of the objects ac- 
creted onto the main branch in different redshift bins. We 
have determined the morphology of our model galaxies by 
using the B-band bulg e-to-disc ratio together w i th the ob- 
servational relation bv lSimien fc de Vaucouleurs! (1 19861 ) be- 
tween this quantity and the galaxy morphological type. For 
the purposes of this analysis, we classify as early types all 



galaxies with AM < 0.4 (AM = M bu i ge - M tota i), as late 
types all galaxies with AM > 1.56, and as lenticulars (SO) 
all galaxies with intermediate value of AM. In Fig. 1101 the 
red indicates early-type morphologies, the green is for SO and 
the blue is for late types. For each redshift bin, the left-most 
histogram refers to normalisation by mass, while the right- 
most histogram refers to normalisation by the number of 
objects. The number of objects accreted in each redshift bin 
is indicated on the top of the histograms. The figure clearly 
shows that, with the exception of the highest redshift bin 
considered, the majority of the objects accreted onto the 
main branch, already have an early-type morphology. 

We saw previously that a BCG is not well described by 
the monolithic approximation. Instead, at z > 0, a BCG 
breaks into an ensemble of distinct galaxies, and it is all of 
them at the same time. A question relevant both to the the- 
oretical and observational sides is then: what is the spatial 
extent of the progenitor set as a function of time? In other 
words, how far do the progenitors of the BCG extend? Fig. 
[TJJshows the spatial distribution of progenitors of five model 
BCGs at different redshifts. Only progenitors more massive 
than 10 10 /i _1 Mq are shown, with symbols that scale with 
mass as in Fig. [JJ The symbols are colour-coded by the star 
formation rate averaged over a time-step (~ 200 Myr). The 
number of progenitors more massive than 10 10 ft -1 Mq is in- 
dicated in the upper-right corner of each panel. The gray 
area marks a 2 x 2 Mpc 2 comoving region and is centred 
on the main progenitor of the BCG. The typical extent of 
the progenitors goes from ~ 10/i _1 Mpc (comoving) at z — 3 
down to about ~ l/i _1 Mpc at z = 0.5. Several panels (e.g. 
middle panels down to z ~ 1) show that the main progenitor 
of the BCG is not necessarily the most massive progenitor 
at that redshift. 

Fig. [12] shows the stellar mass distribution of progeni- 
tors of redshift zero BCGs at redshifts 0.5 (blue), 1 (red), 2 
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Figure 11. Spatial distribution of the progenitors of five model BCGs at redshifts ~ 0.5, 1, 2 and 3. Progenitors more massive than 
10 10 h~ 1 Mq are shown using symbols colour-coded with star formation rate and with size scaling with stellar mass. The number in the 
upper-right corner of each panel indicates the number of such progenitors. The gray area marks a 2 X 2 Mpc 2 comoving region centred 
on the main progenitor of the redshift zero BCG. 
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Figure 12. Stellar mass distribution for the progenitors of local 
BCGs at redshifts 0.5 (blue), 1 (red), 2 (green), and 3 (orange). 
The black histogram shows the stellar mass distribution of local 
BCGs. 



(green), and 3 (orange). The black histogram shows the stel- 
lar mass distribution of local BCGs. The figure clearly shows 
that the mass distribution of progenitors of local BCGs be- 
comes broader and more skewed towards lower masses with 
increasing redshift. At redshift ~ 1 the massive tail of the 



progenitors of local BCGs are about a factor 3 smaller than 
redshift zero BCGs. This factor rises to about 10 at z ~ 3. 



6 THE EVOLUTION OF HIGH-Lx BCGS 

In the previous sections, we discussed one facet of the evolu- 
tion of BCGs, namely: how did low redshift BCGs form? We 
now want to address this issue from an observational per- 
spective. For example, we want to ask if the low assembly 
times measured for our local sample of BCGs are in agree- 
ment with observational results for high-Lx clusters. This 
raises two complementary questions: 

• How do the properties of high-z BCGs differ from those 
of local ones? In other words, what is the observable evolution 
of BCGs ? 

• What is the hierarchical relationship (if any) between 
BCGs observed at different redshifts? In other words, are 
high-z BCGs progenitors of low- z BCGs and are low-z BCGs 
descendants of high- z BCGs? 

We extend the sample used in the previous sec- 
tion by selecting the 125 central galaxies lying in the 
most massive haloes at 6 different redshifts : z ~ 
{0.2; 0.5; 0.75; 1; 1.5; 2}. Our selection corresponds to a 
selection in mass with lower mass limits equal to 
M h /(1O 14 M ) = {5.4; 3.7; 2.7; 1.9; 1.2; 0.67} respectively. 
We chose this somewhat arbitrary mass selection so as to 
preserve the same number density of BCGs at all redshifts. 
This choice would correspond to the (incorrect) monolithic 
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Figure 13. Evolution of the rest-frame if -band absolute magni- 
tude of BCGs as a function of redshift. Black points show model 
BCGs (a small random spread has been added along the x-axis 
for clarity), and the red symbols refer to BCGs of haloes more 
massive than 10 15 Mq. In the upper panel, the dashed lines shows 
predictions from a single burst model with redshift of formation 
5. The lower line is normalised to the median magnitude of lo- 
cal BCGs, and the upper one refers to a normalisation fainter 
by a factor 2. In the lower panel, black points show the resid- 
uals with respect to the model with median normalisation. The 
shaded (green) area shows the rough location of observations (see 
text for details). 



assumption in which each local BCG has a well-defined dom- 
inant progenitor at each redshift. Given that we select the 
most massive haloes at each time-step, our sample is qualita- 
tively comparable to selecting the highest luminosity X-ray 
clusters at each redshift down to a fixed comoving abun- 
dance. 



6.1 Observable evolution 

In the upper panel of Fig. 1131 we show the rest-frame K- 
band absolute magnitude distribution of the (low- and) high- 
z model BCGs. Each galaxy is represented with a small 
random offset with respect to its redshift, for clarity. The 
lower dashed line in this panel shows the magnitude de- 
cline predicted by a single burst model where all the stars 
are formed at z — 5 and evolved passively to the present 
day. This line is normalised to the median magnitude of 
the local sample. The upper line refers to a normalisation 
fainter by a factor 2. Th e model has been built using the 
iBruzual fc Charlotl (|2003T ) code with a solar metallicity stel- 
lar population and a Salpeter IMF with lower and upper 
mass cutoffs mi =0.1 Mq and m u = 100 Mq. In the lower 
panel of Fig. 1131 we show the residuals of the K-band mag- 
nitudes of our model BCGs with respect to the single burst 



model (median normalisation), as a function of redshift. In 
both panels, the red symbols highlight BCGs of haloes more 
massive than 10 15 Mq . Using the local observ ational relation 
between M200 and Lx dPopesso et al.l l2005i). this mass cut 
ro ughly corresponds t o the Lx > 2.3 x 10 44 ergs -1 cut used 
bv lBurke etafl (|2000h . 

As mentioned before, our low-z BCGs cover the same 
range of absolute K-band magnitudes as measured for high- 
Lx clusters. Model data shown in the top panel of Fig. [13] 
indicate that the predicted mass growth between z ~ 1 and 
2 = is of about a factor 3. We note that this is the same 
factor measured in Fig. [7] following the main branch of local 
BCGs. This appears as a coincidence since high- z BCGs are 
a priori different from progenitors of local BCGs (see next 
subsection). This predicted mass growth also seems to be 
in conflict with observational results. The latter are, how- 
ever, based on high-Lx clusters, roughly more massive than 
10 15 Mq. Red symbols in the top panel of Fig. H3l show that 
selecting haloes above a fixed mass cut corresponds to selec- 
tion of "rarer" systems at higher redshifts. For model BCGs 
sitting in these haloes, the mass growth is less than a factor 
2 between z = 1 and z = 0. Although based on low-number 
statistics, this result is compatible with that inferred from 
observational data. 

In the bottom panel o f Fig. 1131 we reproduce the anal- 
ysis bv lBurke et al.l ( |2000| ). The shaded (green) area in this 
panel marks the rough location of their results. The red sym- 
bols again show that a selection above a fixed mass cut gives 
a good agreement with these observational data. The evolu- 
tion of 'rarest' BCGs also appears to be in agreement with 
passive evolution. We find, however, that this does not hold 
if the selection is extended to lower mass systems. 

The above comparis ons should be cons idered only as 
a qualitative indication. iBurke et al.l (2000) use observed 
aperture magnitudes (measured adopting a 50kpc diame- 
ter aperture) and make no attempt to remove the flux due 
to other galaxies falling within the aperture. They also use 
a different model normalisation along with a number of dif- 
ferent details. 

In Fig. [14] we show the rest-frame colour-magnitude 
(AB) relation (upper panel) and the relation between the 
stellar mass and metallicity (lower panel) for our model 
BCGs. Symbols are colour-coded with increasing redshift 
from black to red. The amount of reddening we measure 
between redshift ~ 2 and redshift is ~ 0.2 (considering 
the median colour for each redshift bin), which is very close 
to the amount predicted by the single burst model used 
above. Interestingly the colour-magnitude relation for our 
model BCGs has very little slope at any redshift. The mass- 
metallicity relation shown in the bottom panel of Fig. [14] 
exhibits a more significant, although still weak, slope. 

The results presented in this section change if a differ- 
ent supernovae feedback model is used. Our default model 
in this pape r is th e supernova feedback model used in 
ICroton et al.l (|2006h . in which the ejection rate is propor- 
tional to the star formation rate. We have re peated our anal- 
ysis u sing the supernova feedback model of |Pe Lucia et ail 
(2006), where the ejection rates are computed on the ba- 
sis of energy conservation arguments. This model results 
in less efficient winds in massive galaxies and, as a conse- 
quence in more prolonged star formation activity. While this 
does not significantly change results presented in previous 
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Figure 14. Upper panel : colour-magnitude relation for BCGs as 
a function of redshift. Symbols are colour-coded with increasing 
redshift from black to red. Lower panel : mass-metallicity relation, 
and its evolution as a function of redshift. 



sections, it results in slower evolution of the K-band mag- 
nitudes as a function of redshift which brings more model 
points into agreement with the observational data. However, 
the K-band magnitudes of local BCGs are then predicted to 
be about 0.8 mag brighter and are typically more luminous 
than observed systems. This different feedback scheme also 
produces a larger scatter in the colours and metallicities of 
model BCGs at all redshifts, as well as stellar metallicities 
that are larger than those shown in Fig. [14] by a factor of 
~ 2. The former is a consequence of the more extended 
star formation activity, the latter of the processing of larger 
amounts of gas into stars. 



6.2 Hierarchical evolution 

In a hierarchical scenario, high-z BCGs are not necessarily 
the progenitors of local ones. In this section we want to 
assess the overlap between these two populations. 

A first question is: what is the fraction of high-z BCGs 
which do not end-up in our local sample? In Fig. 1151 we show 
this fraction as a function of redshift. We have highlighted in 
red (resp. blue) the contribution of BCGs residing in haloes 
less (resp. more) massive than the median mass of our sam- 
ples at each redshift, as indicated on the plot. A striking 
remark from Fig. [15] is that ~ 70 per cent of our BCGs at 
z — 2 do not become part of local BCGs, and this does not 
appear to depend strongly on the mass of their host haloes. 
However, this fraction of "lost" BCGs decreases with red- 
shift, as expected, and does so more rapidly for BCGs sitting 
in more massive haloes. We note that high-z BCGs, although 
they do not necessarily end up in our local BCGs, do end 
up in massive galaxies at the centre of massive haloes at 
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Figure 15. Fraction of high-z BCGs that do not end up in our 
sample of local BCGs. The red area shows contribution from 
BCGs residing in haloes with mass lower than the median mass 
measured for the sample at each of the redshifts considered. The 
blue area shows the contribution from massive haloes. 



z = 0. This is illustrated in Fig. 1161 which shows the mass 
of haloes harbouring high-z BCGs versus the mass of the 
haloes harbouring their descendants at z = 0. In this figure, 
open triangles indicate "lost" BCGs, whereas open circles 
indicate BCGs that end-up in our local sample. This fig- 
ure shows that descendents of high-z BCGs simply extend 
towards lower masses our local sample of clusters and we 
emphasise that this does not suggest that "lost" BCGs are 
in any way a different population. 

Turning the question around, we can assess whether the 
massive progenitors of our local sample are different from 
high-z BCGs. Fig. 1171 shows the stellar mass distribution of 
high-z BCGs (red histograms), and that of the progenitors 
of our local BCGs (blue histograms). These distributions are 
both normalised to the number of BCGs at each redshift (i.e. 
125). Overlap of the histograms does not imply that they in- 
clude the same galaxies, but just that objects of given mass 
have equal number densities. At redshift lower than ~ 0.7, 
model BCGs are 'assembled' (see Fig. [6]), and indeed, the 
main branch stands out at the massive end of the blue his- 
tograms. There is then a rather good overlap between the 
mass distributions of high-z BCGs and the massive progen- 
itors of local BCGs. Moving to higher redshift, the massive 
tail of progenitors disappears. As we showed, this is because 
typically at z > 0.7 no branch has yet become dominant in 
the trees of local BCGs. At these redshifts, high-z BCGs still 
have a similar mass distribution as the massive progenitors 
of local BCGs, although Fig. H5l shows they are not the same 
objects. 

At first sight, it may be surprising that our high-z BCGs 
do not end up in our local sample. However, we showed that 
they do belong to the same population as progenitors of local 
BCGs, and do end up as central galaxies of massive haloes. 
The evolution seen in Fig. [K]then appears to be essentially 
driven by the evolution of dark matter haloes. The most 
massive haloes at z = 2 are not in general the progenitors 
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Figure 16. Halo mass distributions. The y-axis shows the masses 
of haloes hosting BCGs at the redshift they are observed; the x- 
axis shows the masses of haloes containing their descendants at 
z = 0. Triangles indicate "lost" BCGs and open circles indicate 
BCGs that end-up in local BCGs. The colour-coding indicates 
redshift, decreasing from purple to red. 



of the most massive haloes at z = and therefore the BCGs 
they contain just fall out of our local sample. 



7 CONCLUSIONS AND DISCUSSION 

In this study, we have used a combination of high-resolution 
iV-body simulations and semi-analytic techniques to inves- 
tigate the formation and evolution of BCGs in a cosmolog- 
ical context. This allows us to explicitly take into account 
the full hierarchy of dark matter structure growth, and its 
impact on galaxy formation. We show that in a CDM Uni- 
verse, this hierarchical build-up implies that central galaxies 
of massive haloes have complex merging histories that are 
not well described by monolithic approximations. 

In the first part of our paper, we have studied in de- 
tail the formation and assembly histories of local BCGs. We 
find that the stars in these systems start forming very early 
(about 50 per cent of the stars in these systems have already 
formed by redshift ~ 5) although the identity of the BCGs 
themselves is not defined before z ~ 0.7. A monolithic ap- 
proximation clearly fails to describe these objects for most 
of their history. Most of the stars are formed in separate en- 
tities in a quiescent mode and most of the mass in the local 
BCGs is gained through accretion of a relatively low num- 
ber (typically 10) of objects more massive thanlO 10 /i _1 M Q . 
Most of these accreted objects, with the exception of those 
accreted at very high redshift, have very low gas fractions 
and star formation rates, quite red colours and an early-type 
morphology. Interestingly, the main branches of BCGs grow 
fast enough that very few major mergers occur along them, 
most of the accretion events being minor mergers. 

It is worth mentioning that our quantitative results do 
depend strongly on the implemented feedback models, both 
from AGNs and from supernovae. T he AGN feedback model 
employed here (jCroton et al.l [20061 ) is extremely efficient in 



switching off cooling in relatively massive haloes even at high 
redshifts. At the s ame time, the super nova feedback model 
we use (also from ICroton et al.l 120061 ) drives strong winds 
that blow out all the gas from galaxies on short timescales. 
The combination of these two processes causes the early 
shutdown of star formation in the progenitors of local BCGs. 
These two particular feedback implementations are among 
the most efficient proposed so far in the literature. The work 
presented here leaves room for less drastic alternatives. 

In the second part of our work, we have studied the ob- 
servable evolution of BCGs by selecting the 125 most mas- 
sive haloes at different redshifts. From these samples of high- 
z BCGs, we measure a mass growth of a factor of ~ 3 from 
z = 1 to z = 0. This is by coincidence very close to the 
growth we find following the main branch of local BCGs. 
These mass growth factors are higher than those inferred 
from observational data. However, when we use a selection 
that more closely matches that of high-Lx clusters, we find 
milder evolution, consistent with the observations. Interest- 
ingly we find that our model BCGs exhibit, at all redshifts, 
a weak correlation between stellar mass and stellar metal- 
licity and an even weaker correlation between their magni- 
tudes and their colours. This changes slightly if a different 
supernova feedback model is employed. Our high-z BCGs 
do end up as central galaxies of massive haloes at z = 
but, because of our cluster sample selection criterion, many 
of their descendants are not included in our local BCG sam- 
ple. However, we have shown that there is no substantive 
physical difference (1) between local BCGs and the descen- 
dents of high-z BCGs, and (2) between high-z BCGs and 
the most massive progenitors of local BCGs. 
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